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Summary 

A  new  third  term  is  incorporated  within  the  multiplicative  decomposition  of  the  deformation 
gradient  in  the  context  of  geometrically  non-linear  mechanics  of  defective  elastic-plastic  crystals. 
This  enhanced  description,  when  applied  to  an  element  of  material  of  finite  volume,  accounts 
explicitly  for  average  local  residual  lattice  distortion  due  to  defects  within  that  volume.  The 
magnitude  of  the  distortion  from  this  third  term,  determined  analytically  for  an  elastic  cylindrical 
volume  of  outer  radius  R  containing  a  single  dislocation  line  threading  its  centre,  is  estimated 
as  [[b/(nR)\^  +f^  -I-  where  b  is  Burgers  vector  magnitude,  /  accounts  for  elastic  non¬ 

linearity,  and  g  accounts  for  core  effects.  For  a  straight  screw  dislocation  in  a  third-order  isotropic 
elastic  medium,  at  a  dislocation  density  of  10  per  cent  of  theoretical  maximum,  bl(jtR)  is  on 
the  order  of  0.1,/  on  the  order  of  0.01  and  g  is  proportional  to  pressure  exerted  by  the  core  and 
can  be  significant.  Predictions  of  stresses  and  dislocation  density  under  simple  shear  and  uniaxial 
compression  demonstrate  differences  from  those  of  usual  crystal  plasticity  at  large  strain  and  for 
high  hardening.  Besides  offering  a  natural  and  precise  delineation  of  contributions  from  dislocation 
velocity  and  dislocation  generation  to  irreversible  deformation,  the  three-term  model  allows  for 
residual  elastic  strains — including  dilatation  observed  in  experiments  and  atomic  simulations — not 
addressed  by  conventional  two-term  crystal  plasticity. 


1.  Introduction 

The  decomposition  of  the  total  deformation  gradient  into  a  product  of  terms  is  a  standard  approach  in 
finite  plasticity  theory.  Such  a  decomposition  into  two  (or  more)  terms  is  re-examined  in  the  present 
work.  Advancements  are  suggested  here  in  an  effort  to  refine  the  physical  meaning  of  each  term  in 
the  decomposition  and  associated  contributions  to  internal  energy  of  the  material. 

Let  F  denote  the  total  local  deformation  gradient  for  an  element  of  crystalline  material.  For 
simplicity,  thermal  effects  are  omitted,  gliding  dislocations  are  the  only  kind  of  defect  considered, 
and  attention  is  restricted  to  a  single  crystal.  Most  typically,  a  two-term  multiplicative  decomposition 
of  the  form 

F  =  (1.1) 

is  invoked,  where  F^  quantifies  elastic  lattice  stretch  and  rotation,  and  F^  accounts  for  plastic  slip 
due  to  dislocation  glide.  This  decomposition  was  perhaps  first  written  explicitly  by  Bilby  et  al.  (1) 
and  was  advanced  soon  after  by  Kroner  (2)  in  a  more  comprehensive  paper.  Several  other  early 
continuum  theories  of  finite  plasticity  with  dislocation  density  were  compared  in  (3). 


^<john.d.claytonl.civ@mail.mil> 


Q.  Jl  Mech.  Appl.  Math,  Vol.  67.  No.  1.  Published  by  Oxford  University  Press  2014. 
This  work  is  written  by  a  US  Government  employee  and  is  in  the  public  domain  in  the  US. 

Advanee  Aecess  publication  31  January  2014.  doi:10.1093/qjmam/hbt026 


128 


J.  D.  CLAYTON 


Three-term  decompositions,  written  as 


F^F^F^F^,  (1.2) 

have  also  been  proposed,  where  the  mapping  F^  accounts  for  physics  distinct  from  recoverable 
elasticity  (F^)  and  lattice  preserving  slip  (F^).  Such  a  decomposition  was  perhaps  first  written 
explicitly  by  Kroner  (2)  and  Bilby  et  al.  (4),  and  was  first  considered  from  a  thermodynamic  point  of 
view  by  Kratochvfl  (5).  In  more  recent  years,  such  a  three-term  decomposition  has  been  introduced 
in  a  number  of  works,  with  various  definitions  or  derivations  offered  for  F^  (6)-(14). 

A  general  premise  (5)  of  three-term  decompositions  in  the  context  of  dislocation  mechanics  is 
that  F^  accounts  for  microstructural  rearrangements  that  increase  the  internal  energy  of  the  material 
element  relative  to  that  of  an  element  that  has  undergone  slip  but  contains  no  dislocations  or  other 
defects  within.  Such  energy  would  be  associated  with  locally  heterogeneous  residual  elastic  strain 
fields  and  atomic  phenomena  (for  example,  core  effects)  not  resolved  by  F^.  Recent  models  for 
metal  plasticity  (6,  10,  11)  have  linked  a  term  akin  to  F^  to  anisotropic  strain  hardening  associated 
with  defects  within  the  element. 

This  article  develops  a  new  non-linear  continuum  description  using  three-term  decomposition 
(1.2),  with  a  focus  on  the  origin  of  each  term  as  it  relates  to  the  underlying  crystal  lattice. 

2.  Conceptual  arguments 

With  reference  to  decomposition  (1.2),  the  following  definitions  are  now  stated  for  a  local  volume 
element  of  crystalline  material  bounded  by  surface  E : 

F  :  deformation  due  to  traction  applied  to  element  boundary  E , 
associated  with  configuration  change  B  B', 

F^  :  deformation  due  to  dislocations  currently  inside  the  element, 
associated  with  configuration  change  B  ^  B\ 

F^  :  deformation  due  to  dislocations  that  have  slipped  through  the  entire  element, 
associated  with  configuration  change  Bq  B. 

Here  ‘deformation’  is  used  in  an  average  sense  for  the  element  since  atomic  positions  within  the 
element  will  not  change  in  an  affine  manner  if  defects  are  present  or  if  traction  along  E  is  not  uniform. 
Local  configurations  of  the  element  are  as  follows: 

Bq  :  unstressed  perfect  lattice  at  some  (initial)  time  ; 

B  :  unstressed  slipped  lattice  containing  no  defects; 

B  :  externally  unstressed  &  internally  stressed  slipped  lattice  containing  defect(s); 

B  :  externally  &  internally  stressed  slipped  lattice  containing  defect(s). 

Representative  configurations  and  deformation  mappings  for  a  cubic  lattice  with  edge  dislocations 
are  illustrated  in  Fig.  1 . 
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Fig.  1  Deformation  mappings  and  configurations  for  a  single  volume  element  within  a  single  crystal 

Regarding  reference  configuration  Bq,  it  is  assumed  that  the  volume  element  is  embedded  within 
a  larger  slab  of  homogeneous  material  that  is  also  unstressed  in  its  (global)  reference  configuration. 
In  this  case,  the  element  and  the  slab  each  consist  of  atoms  arranged  regularly  on  a  perfect  crystal 
lattice  in  the  reference  configuration.  Each  volume  element  fits  exactly  with  its  neighbours  such  that 
slab  in  the  global  reference  configuration  is  simply  connected. 

A  first  local  intermediate  configuration  B  of  the  element  results  from  the  passage  of  dislocations 
completely  through  the  element.  When  separated  from  neighbouring  elements,  this  element  exhibits 
a  perfect  lattice,  but  with  slip  steps  evident  on  E  where  dislocations  have  entered/exited.  The  global 
configuration  corresponding  to  B  is  generally  disconnected. 

A  second  local  intermediate  configuration  B  of  the  element  accounts  for  defects  within  that  element 
that  have  not  yet  fully  passed  through.  This  element  is  free  from  traction  along  its  external  surface  E, 
and  thus  is  in  a  state  of  self  stress.  The  remainder  of  the  slab  may  contain  defects  and  similar  residual 
stresses;  each  element  comprising  the  slab  must,  in  general,  be  mechanically  removed  from  its 
neighbours  so  that  it  too  is  in  a  state  of  self  stress.  External  boundaries  of  the  entire  slab  are  also 
necessarily  traction  free  in  this  global  intermediate  configuration. 

A  local  current  configuration  B  of  the  element  is  achieved  when  the  element  is  stressed  and 
deformed,  along  with  its  neighbours,  so  that  all  volume  elements  fit  together  and  the  body  is 
simply  connected  in  the  global  current  configuration.  In  the  local  current  configuration,  traction 
must  generally  be  applied  to  E .  External  surfaces  of  the  entire  slab  can  support  traction  in  the  global 
current  configuration,  but  need  not  do  so.  If  the  remainder  of  the  slab  outside  of  the  volume  element 
contains  defects,  local  traction  may  not  necessarily  vanish  on  E  in  the  current  configuration  as  a 
result  of  long-range  stress  fields  of  these  defects,  even  if  the  body  is  globally  free  of  applied  surface 
and  body  forces. 

The  preceding  developments  are  applied  to  a  volume  element  containing  a  single  edge  dislocation 
in  Section  3  and  a  single  screw  dislocation  in  Section  4.  Generalisation  to  multiple  defects  follows 
in  Section  5.  Representative  calculations  using  a  constitutive  model  incorporating  the  present 
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three-term  description  are  reported  in  Section  6  for  several  simple  loading  paths,  followed  by  general 
discussion. 


3.  Edge  dislocation 

Consider  a  single  edge  dislocation  in  an  isotropic  elastic  cylinder,  as  shown  in  Fig.  2.  In  an  undeformed 
state  on  the  left,  the  cylinder  is  of  length  L  and  radius  R.  Let  r  denote  the  radius  of  the  image  of 
the  dislocation  core  in  this  undeformed  state.  Coordinates  for  the  left  side  of  Fig.  2  are  {2f ,  Y ,  Z]  — 
{p  cosd,  p  sin0,  Z};  cylindrical  coordinates  are  {p,  6,  Z},  where  6  e  {—it,  it}.  Boundaries  of  the 
outside  of  the  cylinder  and  the  core  are  E  and  a,  and  F"*"  and  F”  are  boundaries  of  the  slip  plane 
along  {Z,  O"*",  Z}  and  {X,  0“,  Z}  for  —R  ^  X  ^  —r.  Assume  L  ^  so  that  the  problem  becomes 
two  dimensional  (plane  strain  with  z  —  Z).  Let 

=  0+)-x(Z,  0“)  =  m(Z,  0+)-m(A,  0“),  -Ri^X^-r.  (3.1) 

The  magnitude  of  the  Burgers  vector  is  b,  and  b  ^  r  <  R.  The  displacement  field 


u  —  x{p,9)  —  X,  v  —  y{p,9)  —  Y. 


(3.2) 


Assume  that  E  and  a  are  traction  free,  and  let  v  denote  Poisson’s  ratio.  The  displacement  field  of 
the  exact  solution  to  this  boundary  value  problem  in  the  context  of  linear  elasticity  is  (15) 


27t 


9 


1 


4(1  -  V) 


(1  -t-  fll  —  02)  sin  29 


fll(p)  = 


8jr(l  —  v) 

(3  —  4y)p^ 
R^  +  r^  ’ 


[2(1  —  2v)lnp  -f  03  -t-  (1  -t-  01  —  02)  cos  20] ; 


aiip)  = 


p2(I?2  ^  ^2)  ’ 


asip)  = 


2p2 


r2  _|_  j.2 


(3.3) 

(3.4) 

(3.5) 


Non-zero  stress  components  (physical  cylindrical  components)  are  Tpp,  Tee,  Tpe  and  = 
v{Tpp  +  Tee)-  When  v  ^  0,  replacement  of  the  plane  strain  boundary  condition  with  free  end 


Fig.  2  Edge  dislocation  parallel  to  Z-axis  centred  in  elastic  cylinder  of  radius  R  and  (out-of-plane)  length  L; 
Burgers  vector  magnitude  is  6;  core  radius  is  r 
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conditions  along  Z  =  -iiL/l  will  lead  to  non-zero  z-displacement  w  and  perturbation  of  the  other 
components  u  and  v. 

Often  it  is  assumed  that  the  cylinder  is  very  large  relative  to  the  radial  coordinate  such  that 
displacements  (and  stresses)  vanish  as  R  ^  oo  giving  ai,  03  — >  0.  Another  frequent  assumption  is 
that  effects  of  the  traction  free  core  can  be  neglected  or  that  the  radial  coordinate  is  sufficiently  larger 
than  the  core  radius,  leading  to  -»  0  and  02  0.  When  such  terms  corresponding  to  outer 

and  core  radii  are  omitted,  strain  energy  per  unit  length  e  is 


e  — 


Gb^ 

471(1  —  v) 


R 

In  —  =  e  +  e  = 
r 


Gb^ 

’l-p-2y2’ 

R  Gb^ 

In _ 1 _ 

2  —  2v  +  2v^ 

\2n 

(1  -  v)2 

111 

r  127r 

(1  -  y)2 

(3.6) 


where  G  is  the  linear  elastic  shear  modulus,  e  is  the  volumetric  contribution  to  the  total  line  energy 
and  e'  is  the  deviatoric  contribution. 

To  motivate  later  discussion  regarding  distortion,  define  the  following  surface  integral: 


ol  def  1  f 


)Ndi:, 


(3.7) 


where  V  —  jtR^L  is  the  volume  of  the  undeformed  domain  (including  the  core  region)  and  u  — 
uei  +  ve2  +  WP3  with  {p,  }  the  Cartesian  basis.  The  unit  outward  normal  to  E  is  A  =  cos  9ei  +  sin  9e2- 
Applying  this  definition  with  u  from  (3.3),  defining  C  =  (1  -|-  ai  —  fl2)lp=R>  and  noting  that  = 
RLd9, 


b  b 

— ^  /  9  sin  9d9  H - _  , 

Ijt^R  8(1  -  v)tx^R  J- 


2jt^R 

b 


(sin^  —  ^  cos  0)1^^  + 


=  — ::  /  u^inOdO 

■7Z 

C  sin  9  sin  29d9 
(3sin6»  -  sin  29)\^^ 


f 


Cb 


48(1  -  v)tt'^R 


,  -{9cos9)\^„  =  — . 
2n^R  jtR 


It  can  be  verified  by  direct  integration  that  other  planar  components  vanish, 


’  —7Z 


(3.8) 

(3.9) 

(3.10) 

(3.11) 


and  since  w  =  0  and  A  ■  ^3  =  0  on  E,  out-of-plane  components  all  vanish  identically:  (/lJ3)jj„  = 
(^23)lin  ~  (^3l)lin  ~  (^Dlin  ~  (^33)lin  ~ 
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If  the  domain  of  integration  is  extended  to  include  the  core  surface  a  and  the  slip  plane  F*,  then 
Vq  —  7t(R^  —  r^)L  is  the  reference  volume  of  a  simply  connected  domain  and 


V 


—  /  u(—sin0)d<7 - 


-  f  udr  +  —  f 
Vq  Jr+  Vo  Jr 


Vo  Jr- 


udr 


bR 


br 


1 


b  b 

jt(R^  —  r^)  7t(R^  —  r^)  jt(R^  —  r^)\_2^  ^  2^ 


b 


n{R  +  r)  n{R  +  r) 


=  0, 


(3.12) 


consistent  with  the  proof  that  in  the  context  of  linear  elastostatics,  volume-averaged  strain  ^ 

N  +  N  ®  u)dS  —  2^  f[S/u  +  ( VM)^](iVo  must  vanish  in  a  domain  with  uniform  material  properties 
that  is  in  a  state  of  self  stress  (that  is,  a  homogeneous  body  with  no  traction  applied  to  external 
surfaces,  traction  continuity  across  internal  surfaces,  and  no  applied  body  forces)  (16).  Summed 
contributions  from  core  and  slip  surfaces  (ct  and  F^)  to  expressions  analogous  to  (3.12)  for  other 
planar  components  of  distortion  also  vanish. 

Three-term  decomposition  F  —  F^F^F^  of  (1.2)  is  now  revisted  in  the  context  of  a  volume 
element  containing  an  isolated  edge  dislocation  embedded  in  a  larger  slab  of  material  that  may 
contain  other  defects  in  its  deformed  state.  Consider  the  cylindrical  volume  shown  in  Fig.  2,  which 
is  identihed  with  a  volume  element  of  material  to  which  (1.2)  is  applied. 

The  first  local  intermediate  configuration  B  of  this  element  is  identified  with  the  left  part  of  Fig.  2: 
no  local  stress  within  the  element  or  traction  along  its  boundaries.  The  elastic  cylinder,  minus  the 
core  region,  is  simply  connected  since  surfaces  F^  are  treated  as  distinct. 

The  second  local  intermediate  conhguration  B  of  the  element  is  identihed  with  the  right  part  of 
Fig.  2.  This  element  is  free  from  traction  along  its  external  surface  E,  and  thus  is  in  a  state  of  self 
stress  in  two  dimensions  (recall  out-of  plane  stress  Tzz  is  associated  with  the  Poisson  effect  and  plane 
strain).  The  elastic  cylinder,  when  the  core  region  is  removed,  now  becomes  multiply  connected  when 
referential  surfaces  F^  are  displaced  relative  to  one  another  by  the  Burgers  vector  and  then  bonded 
together.  The  remainder  of  the  slab  may  contain  defects  and  similar  residual  stresses;  each  element 
comprising  the  slab  must,  in  general,  be  removed  from  its  neighbours  so  that  it  too  is  in  a  state  of 
self  stress. 

A  local  current  conhguration  B  of  the  element  is  achieved  when  the  element  is  stressed  and 
deformed,  along  with  its  neighbours,  so  that  all  volume  elements  ht  together  and  the  body  is 
continuous  in  the  global  current  conhguration.  In  the  local  current  conhguration  for  each  element 
(not  shown  in  Fig.  2),  traction  is  generally  non-zero  on  E. 

First  consider  plastic  deformation  map  F^,  which  accounts  for  n  dislocations  that  have  already 
moved  completely  through  the  cylinder,  from  left  to  right.  Attention  is  restricted  to  edge  dislocations 
on  the  single  slip  plane  F,  which  is  now  extended  to  bisect  the  entire  cylinder  along  F  =  0.  In  this 
case,  plastic  deformation  is  the  simple  shear 


F^  djf  j  ^ 


(3.13) 


where  the  cumulative  plastic  shear  from  relative  sliding  of  the  two  sections  of  the  cylinder  is 


def 

y  = 


bLR&mOdO  — 


2bn 

jtR 


(3.14) 
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There  is  no  plastic  volume  change,  —  detF^  =  1  +  ej  ■  ^2  =  1-  Since  =  1  —  yei  ®  C2, 
it  follows  that  —  F^F^~^  —  yei  ®  P2-  Assume  for  simplicity  that  the  dislocations  are  equally 
spaced;  their  density  is  then  ^  =  1  / Assume  also  that  dislocations  move  with  constant  velocity 
u  =  2Rh  in  the  direction  of  ej,  where  h  is  the  number  of  dislocations  per  unit  time  that  transit  the 
diameter  of  the  cylinder;  then  the  usual  so-called  Orowan  equation  is  recovered: 


y 


2bh 

jtR 


bv 

jxR^ 


—  l^bv. 


(3.15) 


Next  consider  the  mapping  F^.  In  Cartesian  coordinates,  this  can  be  written 

F^  = 


(3.16) 


The  distortion,  by  definition,  can  be  further  partitioned  into  the  sum 

^'  =  ^Hn+^LHn+^Le-  (3-17) 

The  linear  term  is  defined  as  the  distortion  of  the  outer  boundary  E  of  the  elastic  cylinder  calculated 
from  linear  elasticity,  where  from  (3.7)  and  (3.8), 

^lin  =  ^  Cl  ®e2-  (3.18) 

This  non-symmetric  term  clearly  includes  both  stretch  and  rotation.  The  other  two  terms  are 
corrections  accounting  for  non-linear  elastic  effects  and  core  effects  that  are  omitted  by  the  linear 
elastic  solution.  Because  of  non-linearity,  these  terms  in  general  cannot  be  determined  exactly  from 
superposition  of  solutions  of  decoupled  boundary  value  problems. 

The  non-linear  correction  is  defined  as  (17) 


(3-19) 

or  in  indicial  notation 

where  Ua  are  local  displacements  of  the  self-equilibrated  body  in  configuration  B  measured  fromZ^, 
which  are  local  coordinates  in  configuration  B.  Let  Ca^yS  denote  the  second-order  elastic  constants 
and  Sa^ys  the  complementary  elastic  compliance.  Let  CapyS^t;  denote  third-order  elastic  constants, 
omitted  in  the  preceding  linear  elastic  analysis  of  the  edge  dislocation.  Particular  forms  of  elastic 
stiffness  tensors  and  representative  values  for  isotropic  solids  are  described  in  Appendix  A.  Self¬ 
equilibrium  conditions  for  the  cylinder  in  configuration  B  lead  to  vanishing  of  its  volume  averaged 
Cauchy  stress,  giving  (16, 17) 


(^nonlin)^y  I 


r  dup  duy  f  dua  9m; 

Jv  dX,  dXs  Jv  dX,  dx. 

r\  I  /»  A 


filJ^  filJ^ 


(3.21) 


to  second  order  in  displacement  gradients  Vh,  assuming  elastic  moduli  are  independent  of  location 
within  V.  Terms  of  orders  three  and  higher  in  displacement  gradients,  involving  elastic  constants 
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of  orders  four  and  higher,  are  omitted.  By  this  dehnition,  is  symmetric  (that  is,  due  to  the 

symmetry  of  the  first  two  components  of  the  elastic  compliance  tensor);  a  rigid  body  rotation  could  be 
superposed  without  violating  the  self-equilibrium  constraint.  Logically,  vanishes  in  the  linear 

elastic  approximation  since  products  of  order  two  in  displacement  gradients  are  inconsequential  in 
the  stress  field  of  linear  elasticity. 

Evaluation  of  integral  (3.21)  requires  knowledge  of  the  displacement  field  u{X)  for  the  non-linear 
elastic  solution  of  an  edge  dislocation  embedded  in  a  cylinder  of  a  third-order  elastic  material  with 
finite  boundaries.  No  known  closed-form  solution  exists  to  this  boundary  value  problem  that  permits 
ready  evaluation  of  the  integral.  Zubov  (18)  tersely  presents  an  exact  analytical  solution  for  an  edge 
dislocation  embedded  in  an  infinitely  extended  domain  of  a  semi-linear  material;  the  displacement 
field  of  this  solution  is  expressed  in  terms  of  complex  functionals  that  cannot  be  easily  inserted  into 
(3.21).  The  stress  field  of  the  non-linear  solution  does,  however,  demonstrate  significant  differences 
from  the  linear  elastic  solution  at  distances  up  to  the  magnitude  of  several  Burgers  vectors  from  the 
core.  Seeger  et  al.  (19)  give  an  approximate  iterative  solution  for  an  edge  dislocation  accounting  for 
up  to  second-order  terms  in  displacement  gradients,  but  the  displacement  field  is  written  as  a  series  of 
complex  functions  that  do  not  permit  analytical  evaluation  of  (3.21).  On  the  other  hand,  a  tractable, 
yet  still  approximate,  analytical  non-linear  elastic  solution  does  exist  for  the  screw  dislocation,  as 
will  be  discussed  later  in  Section  4. 

An  analytical  approximation  for  the  dilatation  due  to  non-linear  elastic  effects  from  residual 
stresses  does  exist,  however,  for  isotropic  materials  as  well  as  several  cubic  crystal  classes.  When 
contributions  of  products  of  orders  three  and  higher  in  local  displacement  gradients  are  omitted, 
the  volume  change  of  an  isotropic  cylinder  with  an  edge  dislocation  due  to  non-linear  elasticity 
is  (16, 17) 
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h^^ln- =Ah2^1n-,  (3.22) 

r  r 


where  K  is  the  linear  elastic  bulk  modulus,  K'  and  G'  are  pressure  derivatives  of  the  incremental  bulk 
and  shear  moduli  in  the  stress-free  reference  state  (and  can  be  related  to  third-order  elastic  constants 
as  described  in  Appendix  A),  A  is  a  constant  dimensionless  function  of  the  elastic  constants,  and 
^  =  L/V.  In  this  solution,  e  correctly  incorporates  only  contributions  from  second-order  elastic 
constants  (16),  but  effects  of  non-linearity  and  boundary  conditions  on  the  dislocation  core  and  free 
outer  radius  to  this  contribution  to  elastic  energy  are  neglected.  To  integral  terms  of  at  least  hrst 
order  in  displacement  gradients  (17), 

AV/V  ^  det(l  +  -  1  ^  tr^l„„,„.  (3.23) 

Knowledge  of  dilatation  alone  is  insufficient  to  determine  all  components  of  resulting  from 

an  edge  dislocation  in  an  isotropic  cylinder.  Following  (19),  (15),  assume  that  the  displacement  field 


(3.24) 
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where  is  the  linear  elastic  solution  and  is  a  correction  that  accounts  for,  in  the  resulting 
stress  field,  non-linear  elastic  effects  of  up  to  second  order  in  displacement  gradients.  Let  = 
Assume  again  that  the  cylinder  is  long  (plane  strain)  and  that  the  second- 
order  correction  is  independent  of  the  orientation  of  slip  plane  F .  These  symmetry  arguments  suggest 

—  bxip)  cos 0,  —  bx(p)sin0,  —  0,  (3.25) 

where  x(p)  is  a  dimensionless  correction  to  the  radial  displacement  field  from  non-linear  effects  and 
p  is  the  radial  coordinate  of  the  undeformed  cylinder  (left  side  of  Fig.  2).  Applying  the  divergence 
theorem. 


T  1  f  3m(2)  1  f 

(4onHn)ll  -^dV  -yj^bx  COS 

=  \x(R)  r  cos2  Ode  =  ^x(R),  (3.26) 

ttR  R 

T  1  r  3vi2)  1  f  , 

(/ino„lin)22  ^  y  -j^dV  =  y  J^bx  sin 

b  7  b 

=  —xiR)  /  =  -x(R),  (3.27) 

7tR  R 

and  all  other  terms  vanish  identically: 


Therefore,  from  (3.23),  a  coarse  approximation  for  the  edge  dislocation  is 
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(3.29) 


Taking  representative  values  from  Table  A.  1  of  Appendix  A,  (3.22)  and  (3.29)  give  A  ^  and,  with 

and  r  ^  b,  (^J„„ii„)ii  =  ^ 

No  analytical  solution  exists  for  ^l^re-  ^^y  definition,  continuum  elasticity  theory  breaks  down 
inside  the  core  region,  but  atomic  theory  can  be  used  to  estimate  its  contribution.  Term 
accounts  for  uncertainty  in  core  size  and  shape  (20)  and  non-zero  traction  imparted  by  the  core 
on  the  surrounding  continuum  (15).  In  elastic  calculations  considered  thus  far,  assumptions  of  a 
cylindrical  core  with  minimum  radius  r  ^  b  (20)  imparting  no  traction  on  the  surrounding  material 
have  been  used.  Presumably,  for  a  small  cylinder  (large  defect  density),  components  of  could 
be  comparable  or  larger  than  those  of 

In  atomic  simulations  of  an  edge  dislocation  in  a-iron  (21),  total  dilatation  was  found  to  obey 
AF/y  =  [0.37  +  0.07  ln(R/b)]b^^,  with  contributions  from  non-linear  and  core  effects  estimated 
as  similar  orders  of  magnitude,  though  several  arbitrary  choices  were  made  in  partitioning  non¬ 
linear  and  core  contributions.  In  atomic  simulation  of  the  strain  field  of  this  same  dislocation  type 
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(22),  it  was  found  that  deviations  from  the  linear  elastic  Volterra  solution  could  be  fit  to  an  elliptical 
expansion  field  centred  slightly  below  the  centre  of  the  extra  half-plane.  The  procedure  for  extracting 
residual  strain  fields  induced  by  dislocation  lines  in  atomic  simulations  is  intricate,  requiring  pages 
of  discussion  in  (21)  too  lengthy  to  repeat  here  and  to  which  the  reader  is  referred  for  details.  Only 
their  preferred  end  result  (p.  3896  of  (21)  giving  the  just-mentioned  form  of  AV/V)  is  used  in  later 
Table  1;  this  result  may  of  course  be  subject  to  possible  inaccuracies  associated  with  the  numerical 
method  of  relaxation,  imposition  of  non-unique  initial  conditions  and  boundary  constraints,  selection 
of  Eulerian  versus  Lagrangian  coordinates,  choice  of  atomic  potential,  etc.,  noted  in  (21). 

Three  other  obvious  idealisations  have  been  made  in  the  definition  of  F^.  First,  the  linear  term 
would  be  slightly  different  if  the  top  and  bottom  faces  of  the  cylinder  were  relaxed  and  Poisson’s  ratio 
were  non-zero.  Second,  the  dislocation  is  assumed  centred  in  the  cylinder.  More  general  dislocation 
arrangements  are  mentioned  in  Section  5.  Third,  isotropic  elasticity  has  been  assumed.  Anisotropic 
solutions  would  be  more  realistic,  but  general  analytical  solutions  are  not  available  in  compact 
equations  that  are  easily  manipulated,  except  for  cubic  crystals  of  certain  classes  for  which  a  solution 
for  volume  change  similar  to  (3.22)  exists  (16, 17). 

Combining  (3.16),  (3.17),  (3.18),  and  (3.29), 


-  1  +  ^lin  +  ^nonlin  +  ^!:ore 

b  A  R  r 

^  1  +  — ei  ®  ^2  -t-  —  ^  In  -(ei  ®  ei  -f  e2  ®  C2)  +  Pcore 

=  1-1-  {b^^l7r)^^^ei  ®  C2  +  <8>  +  ^2  <S>  C2)  +  (3.30) 


If  the  core  term  is  omitted  and  the  approximations  r  —  b  and  A  =  ^  are  used, 

u2 

—  F^  (b^^  ^62+  ln[(7rh^O~^^^](^l  ®  +  ^2  <S>  ^2)-  (3.31) 

Non-vanishing  components  of  (3.31)  are  shown  in  Fig.  3,  along  with  magnitude 

Ill'll  =  =  mL)nf  +  2[(^Lniin)n]'  + 

^  {[(b/(jrR)f  +  2^[h/(7r/?)]4[ln(/?/h)]2}l/2,  (3.32) 

where  g  is  0(||)3gQjg||)  and  is  again  omitted  in  the  final  approximation.  The  linear  contribution 
dominates.  The  theoretical  maximum  dislocation  density  is  approached  as  R  ^  f^min  =  ^ 
b^^  b^^max  =  ^ i  fo'"  example,  ^max  ^  3  X  10'*m“^  for  ^  0. 1  nm^.  Several  values  are  listed 
in  Table  1 .  Both  linear  and  dilatational  contributions  can  be  significant  relative  to  applied  elastic 
strain  associated  with  F^,  since  elastic  shear  strains  at  yielding  in  most  crystals  are  on  the  order  of 
10“^  —  10“^.  Also  shown  in  the  rightmost  column  of  Table  1  is  the  normal  strain  (estimated  as  half 
of  the  total  dilatation)  reported  from  atomic  simulation  of  a  [100]  edge  dislocation  a-iron  (21).  For 
R/b  ^  1.7,  results  of  atomic  calculation — which  account  for  elastic  anisotropy,  elastic  non-linearity 
and  core  effects — are  no  more  than  one  order  of  magnitude  larger  than  the  third  column  of  Table  1 
that  addresses  isotropic  non-linear  elasticity  alone,  though  generic  properties  from  Table  A.  1  used 
in  calculation  of  the  latter  are  not  intentionally  representative  of  iron. 
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Fig.  3  Analytical  estimate  of  shear  component  /Sj2  and  dilatation  A  V'/(2F)  for  single  edge  dislocation  centred 
in  isotropic  elastic  cylinder  of  radius  R.  Normalized  dislocation  density  is  f /fmax  =  iP'  jR^ 


Table  1  Estimated  contributions  to  for  isolated  edge  dislocation 


R/b 

^  /(imax 

(Ain)l2 

AV/(2V) 

Atomic  simulation  (21) 

1 

1 

0.32 

0 

5.9  X  10-2 

1.7 

0.35 

0.19 

2.4  X  10-3 

2.3  X  10-2 

3 

0.11 

0.11 

1.6  X  10-3 

7.8  X  10-3 

10 

1.0  X  10-2 

3.2  X  10-2 

3.0  X  10-‘* 

8.5  X  10-“^ 

100 

1.0  X  10-"^ 

3.2  X  10-3 

6.1  X  10-® 

1.1  X  10-3 

Effects  of  elastic  anisotropy  are  considered  in  Appendix  C.  There  it  is  shown  that  AV/(2V) 
associated  with  may  be  expected,  due  to  elastic  anisotropy,  to  increase  or  decrease  relative  to 

the  isotropic  solution  by  no  more  than  a  factor  of  around  six  for  dislocations  on  typical  glide  systems 
in  many  Face  Centred  Cubic  (FCC)  or  Body  Centred  Cubic  (BCC)  metals.  Note  that  which 
tends  to  dominate  the  magnitude  of  total  residual  lattice  distortion,  is  unaffected  by  anisotropy. 

From  the  polar  decomposition  theorem  applied  to  (1.2),  F  —  V^R^F^F^ .  Inverse  stretch  is 
associated  with  unloading  of  the  cylinder  from  its  current  deformed  state  (23),  and  is  elastic 
rotation  that  encompasses  all  rigid  body  rotation  in  the  absence  of  defects.  Presuming  F^  and 
F^  are  calculable  from  dehnitions  given  above,  and  that  total  deformation  F  is  imposed,  then 
F^  —  FF^~^F^~^  is  simply  the  remainder. 

Several  points  of  clarihcation  are  in  order.  Yielding  is  indicated  by  non-negligible  plastic 
deformation  described  by  F^,  which  in  turn  depends  on  plastic  slip  y  in  (3.13).  For  monotonic 
loading,  =  1  and  y  =  0  at  stresses  below  the  yield  limit.  Residual  lattice  distortion  can  be  non¬ 
zero  whenever  the  body  contains  dislocations,  regardless  of  whether  or  not  plastic  deformation  takes 
place  concurrently  (that  is,  regardless  of  whether  or  not  dislocations  are  mobile).  It  is  emphasised 
that  F^  is  not  ‘plastic  deformation’  but  rather  a  kind  of  ‘residual  elastic  lattice  deformation’,  which 
justifies  the  summation  of  linear,  non-linear  and  core  elastic  helds  in  (3.17)  and  (3.30).  The  solution 
for  the  linear  elastic  fields,  though  well  known,  has  been  included  here  in  order  to  enable  derivation 
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of  in  a  form  not  presented  elsewhere.  None  of  the  components  of  (linear,  non-linear  or  core) 
can  be  included  in  since  by  dehnitions  given  in  Section  2,  F^  consists  of  reversible  or  recoverable 
elastic  deformation,  while  F^  consists  of  residual  deformation.  Instead  of  the  summation  approach 
in  (3.17),  a  multiplicative  split  such  as  F^  =  ^lin^nonlin^core  could  be  assumed.  The  difference 
between  the  two  approaches  would  be  0(  ||)S^  |P),  which  a  posteriori  should  be  small  relative  to  ||)S^  || 
according  to  Table  1.  Therefore,  use  of  such  a  multiplicative  split  is  not  warranted  at  present. 


4.  Screw  dislocation 

Now  consider  an  isolated  screw  dislocation  of  Burgers  magnitude  b  aligned  along  the  Z-axis, 
contained  within  a  cylinder  of  the  same  geometry  as  shown  on  the  left  side  of  Fig.  2.  Boundaries  of 
the  outside  of  the  cylinder  and  the  core  are  S  and  ct,  and  r"*"  and  F"  are  boundaries  of  the  slip  plane 
along  {X,  O"*",  Z}  and  {X,  0“,  Z}  for  —R  ^  X  ^  —r.  Let 

b  —  z{p,Tt)  —  z{p,—Tt),  R^p^r.  (4.1) 


Let  the  resultant  force  and  torque  on  each  end  of  the  cylinder  at  Z  =  ±Z./2  vanish.  In  physical 
components,  the  displacement  field  for  the  linear  solution  is  (15) 


U0  = 


bpZ 
71 R^  ’ 


Uy  —  W  — 


be 

lit 


(4.2) 


The  analogue  of  (3.6),  for  L  — >•  oo,  is 
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Applying  definition  (3.7)  for  this  isolated  screw  dislocation. 
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All  other  components  vanish; 

= ("y,,. ('‘5' 

Three-term  decomposition  F  =  of  (1.2)  is  now  revisted  in  the  context  of  a  volume 

element  containing  an  isolated  screw  dislocation  embedded  in  a  larger  slab  of  material  that  may 
contain  other  defects  in  its  deformed  state. 

First  consider  plastic  deformation  map  F^,  which  accounts  for  n  dislocations  that  have  already 
moved  completely  through  the  cylinder.  Attention  is  restricted  to  screw  dislocations  on  the  single 
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slip  plane  F,  which  is  now  extended  to  bisect  the  entire  cylinder  along  its  axis.  In  this  case,  plastic 
deformation  is  the  simple  shear 

(M  j  ^ 

where  the  cumulative  plastic  shear  from  relative  sliding  of  the  two  sections  of  the  cylinder  is 


def 

y  = 


bL  R  sin  0d9 
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(4.7) 


There  is  no  plastic  volume  change:  —  det  —  I  +  ej,  ■  62  =  1 .  Assume  for  simplicity  that  the 

dislocations  are  equally  spaced;  their  density  is  then  C,  —  Assume  also  that  dislocations  move 
with  constant  velocity  v  —  2Rh\  then 


Y 


2bh 

JtR 


bv 

icR^ 


—  ^bv. 
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Equations  (4.7)  and  (4.8)  are  identical  to  (3.14)  and  (3.15),  respectively. 

Next  consider  the  mapping  F^.  Equations  (3.16)  and  (3.17)  still  apply.  The  linear  contribution  is 


JtR  ^ 


)  62 


(bh/n) 


1/2 


^3  0  ^2  5 


(4.9) 


which  includes  both  strech  and  rotation.  The  non-linear  correction  is  defined  as  in  (3.19),  (3.20)  and 
(3.21).  When  contributions  of  products  of  orders  three  and  higher  in  local  displacement  gradients 
are  omitted,  the  volume  change  of  an  isotropic  cylinder  with  a  screw  dislocation  due  to  non-linear 
elastic  effects  is  (16,  17) 


AT 


9  /?  r.  R 

r^ln-  =Afe2^1n-, 


(4.10) 


where  A  is  a  different  constant  dimensionless  function  of  the  elastic  constants,  and  ^  —  L/V.  Taking 
values  from  Table  A.  1  as  representative,  A  Rs  In  this  solution,  e  —  e’  correctly  incorporates  only 
contributions  from  second-order  elastic  constants  (16),  but  effects  of  nonlinearity  and  boundary 
conditions  on  the  dislocation  core  and  free  outer  radius  to  this  contribution  to  elastic  energy  are 
neglected. 

Knowledge  of  the  dilatation  alone  is  insufficient  to  determine  all  components  of  resulting 

from  a  screw  dislocation  in  an  isotropic  cylinder.  Evaluation  of  integral  (3.21)  requires  knowledge  of 
the  displacement  field  u{X)  for  the  non-linear  elastic  solution  of  a  screw  dislocation  embedded  in  a 
cylinder  of  a  third-order  elastic  material  with  finite  boundaries.  No  known  exact  closed-form  solution 
exists  to  this  boundary  value  problem  that  permits  ready  evaluation  of  the  integral.  Zubov  (24) 
presents  an  exact  analytical  solution  for  a  screw  dislocation  embedded  in  an  incompressible  non¬ 
linear  elastic  material.  The  stress  field  of  this  non-linear  solution  demonstrates  significant  differences 
from  the  linear  elastic  solution  as  p  — )►  7>.  The  incompressibility  constraint  (v  —>■  0.5)  is  not 
realistic  for  most  crystals  and  precludes  calculation  of  dilatation  due  to  geometric  and  material 
nonlinearity.  Teodosiu  (15)  gives  an  approximate  iterative  solution  for  a  screw  dislocation  embedded 
in  a  cylindrical  domain  accounting  for  up  to  second-order  terms  in  displacement  gradients,  where 
the  cylinder  is  infinitely  long  but  may  have  free  boundaries  at  r  and  R.  This  solution  will  be  used  to 
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inform  (3.21)  in  what  follows.  Earlier,  Willis  (25)  obtained  a  similar,  but  not  identical,  approximate 
solution  for  a  screw  dislocation  in  an  infinite  medium  using  Green’s  functions. 

Following  (15)  (p.  211),  assume  that  (3.24)  holds,  where  is  the  linear  elastic  solution  and 
is  a  correction  that  accounts  for,  in  the  resulting  stress  field,  non-linear  elastic  effects  of  up  to  second 
order  in  displacement  gradients.  Let  uP''^  =  The  approximate  non-linear 

solution  is  of  the  form 

—  ri(p)  cos  6,  —  rjip)  sin  6 ,  =  0,  (4.11) 


where  r]{p)  is  a  correction  to  the  radial  displacement  field  from  non-linear  effects,  and  p  is  the  radial 
coordinate  of  the  dislocated  cylinder.  When  the  core  surface  at  p  =  r  is  traction-free,  the  solution  is 
(15)  (p.  231) 

Cl  In  p 

p(p)= - |-C2p-t-C3 - ,  (4.12) 

P  P 

where  c\,  C2  and  C3  depend  on  b,  r,  R,  and  the  second-  and  third-order  elastic  constants,  and  are 
given  in  Appendix  B.  Applying  the  divergence  theorem, 

=  ^  r  -  >imd0,  (4.13) 

7t[R- r](R)Y 


where  the  Lagrangian  radial  coordinate  of  E  is  /?  —  r](R),  since  p  is  the  Eulerian  radial  coordinate 
in  the  solution  of  (15).  Non-vanishing  components  of  (4.13)  are 


^^nonlin^ll 


(^Lnlin)22 


7T[R 


7T[R 


1  n 

-  rim 

1  r 

-  rim  J-n 


9  n(/?) 

ri(R)  cos^  9d6  = 


ri(R)  sin^  9d9  — 


R-riiR)' 

ri(R) 

R-p(Ry 


all  other  terms  vanish  identically,  that  is,  (3.28)  applies  for  the  screw  dislocation. 
Combining  (3.16),  (3.17),  (4.4),  (4.14)  and  (4.15), 


(4.14) 

(4.15) 


-1+  +  /Jjjoniin  +  /3core 

b  p(R)  T 

1  +  —es  ®  C2  +  ®  gl  +  g2  <S>  g2)  +  Pcore 

=  1-1-  Qp-t,  Inf  1^62,  ®  g2  +  ,  1  0  Cl  +  g2  <S>  ^2)  +  (4-16) 


Calculable  terms  of  (4.16)  and  dilatation  Ay/(2y)  from  (4.10)  are  compared  in  Eig.  4,  along  with 
the  magnitude 

11^^  II  =  =  {[(Ain)32]^  +f  + 

^{[(h/(7r/?)]2  +  2[(/i;„„ii„)n]^}'/^ 


(4.17) 
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Fig.  4  Analytical  estimate  of  shear  component  non-linear  contribution  /S  j  j  =  (622  ’  dilatation  A  V/(2  V) 

for  single  screw  dislocation  centred  in  isotropic  elastic  cylinder  of  radius  R.  Normalised  dislocation  density  is 
C/?max  =  fR?'.  In  part  (a),  total  and  linear  curves  are  visibly  indistinguishable. 

Table  2  Estimated  contributions  to  for  isolated  screw  dislocation. 


R/b 

^  /  ^max 

(^Lh2 

^^nonlin^^l 

^V/{2V) 

(^™re)ll  •  G/p 

1 

1 

0.32 

5.9  X  10-2 
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00 

1.7 

0.35 

0.19 

2.7  X  10-2 

9.1  X  10-"^ 

0.40 

3 

0.11 

0.11 

8.2  X  10-3 

6.1  X  10-“^ 

9.3  X  10-2 

10 

1.0  X  10-2 

3.2  X  10-2 

7.0  X  10-"* 

1.2  X  10-"^ 

7.6  X  10-3 

100 

1.0  X  lO-"* 

3.2  X  10-3 

6.2  X  10-® 

2.3  X  10-® 

7.5  X  10-3 

where/  =  =  ^(^nonlin)22  is  the  non-linear  elastic  contribution,  and  g  is  Odl/Sj-o^e  II) 

and  is  omitted  in  the  final  approximation.  In  the  context  of  (4.19)  and  (4.20)  to  be  discussed  later, 
g  —  \/2(/core)ll  =  's/2(/^Qj.g)22  oc  p/G  and  increases  with  increasing  f,  with  p  radial  pressure 
exterted  by  the  core  on  the  surrounding  crystal.  The  linear  term  exceeds  /,  though  the  latter  is  not 
insignificant.  Several  values  are  listed  in  Table  2.  Linear,  non-linear  and  core  contributions  can  all 
be  substantial  relative  to  applied  elastic  strain  associated  with  F^. 

Note  that  an  approximation  of  the  non-linear  contribution  in  terms  of  dilatation  such  as  that  made 
in  (3.29)  is  not  necessary  for  the  screw  dislocation  since  the  second-order  non-linear  elastic  solution 
for  r]  is  available  for  this  geometry.  As  is  clear  from  Table  2,  non-linear  term  (/Jonij^)!!  exceeds 
dilatational  term  A  y/(2y )  by  several  orders  of  magnitude  near  the  theoretical  maximum  density,  but 
both  terms  are  of  comparable  magnitude  for  ^  / (I’max  ^  0.01.  The  former  is  presumably  more  accurate 
because  it  accounts  consistently  for  non-linear  contributions  to  the  elastic  energy  of  the  dislocation 
and  the  free  boundary  at  p  =  /?;  however,  it  is  still  an  approximation  because  non-linear  contributions 
to  the  average  axial  stress  need  not  vanish  (the  cylinder  is  not  self-equilibrated  in  the  axial  direction). 
The  latter  approximation  (dilatation)  relies  on  elastic  energy  per  unit  length  of  a  screw  dislocation 
in  an  infinite  linear  elastic  medium.  Both  approximations  [(/noniin)ll  and  Ay/(2y)]  omit  higher- 
order  corrections  involving  contributions  of  terms  of  order  three  and  higher  in  local  displacement 
gradients  to  stresses  and  average  strains,  and  both  also  neglect  contributions  of  the  dislocation  core  to 
traction  at  p  =  r.  For  values  of  material  constants  considered  here,  both  approximations  do  however 
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demonstrate  positive  volume  change  that  has  been  confirmed  experimentally  for  cold- worked  metals 
with  cubic  symmetry  (26). 

Now  consider  for  edge  dislocation,  continuum  elasticity  theory  cannot  be  applied  to 

describe  the  core  of  a  screw  dislocation,  and  accounts  for  uncertainty  in  core  size  and  shape 
and  non-zero  traction  imparted  by  the  core  on  the  surrounding  crystal.  In  an  isotropic  linear  elastic 
solid,  the  contribution  of  a  radial  pressure  p  acting  on  core  surface  p  =  r  to  the  dilatation  of  the  solid 
can  be  determined  analytically.  Displacement  u  imparted  by  such  a  core  pressure  is  radial  (15): 


m(p)  = 


pr 


2G(R^  -  r2) 
Non-zero  components  of  follow  as 
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(4.18) 


(4.19) 

(4.20) 


Predicted  components  of  can  be  comparable  or  larger  than  those  of  and  the  core 

contribution  becomes  singular  as  the  cylinder  shrinks  to  the  core  itself  (R/b  1).  Representative 
values  are  listed  in  Table  2  for  r  =  b  and  v  =  4 .  Core  pressure  p  at  a  prescribed  radius  r  must 
be  supplied  by  atomic  theory;  recent  calculations  for  dislocations  in  aluminum  (27)  suggest  that 
Ipl  <  O.IG.  Atp/G  =  ^  for  R  ^  1.71?.  At  very  large  pressures  (for  example,  p 

approaching  G)  the  linear  elastic  estimate  of  in  (4.19)  and  (4.20)  would  not  be  expected  to 
be  accurate,  and  non-zero  core  traction  would  also  alter  the  solution  in  Appendix  B  and  Table  2 
for  the  non-linear  elastic  contribution  to  dilatation.  Density  functional  theory  calculations  (28)  for 
core  contributions  to  displacements  for  various  [111]  screw  dislocation  dipole  arrangements  in  a- 
iron  predict  dilatation  perpendicular  to  the  dislocation  line  and  contraction,  smaller  in  magnitude, 
parallel  to  the  line.  These  calculations  incorporate  relaxed  periodic  boundary  conditions,  while  the 
elastic  solutions  used  herein  for  screw  dislocations  omit  the  possibility  of  deformation  parallel  to 
the  infinitely  extended  dislocation  line. 

The  same  comments  regarding  the  polar  decomposition  of  near  the  end  of  Section  3  apply  for 
the  isolated  screw  dislocation. 


5.  Generalisations 


5.1  Multiple  and  moving  dislocations 

Since  and  F^  are  not  state  functions  but  instead  depend  on  load  history  for  arbitrary  deformation 
paths,  it  is  most  appropriate  to  consider  rate  forms  of  these  quantities  for  generalisation  to  an  element 
of  a  single  crystal  with  multiple  dislocations  whose  number  and  position  change  with  time.  Let  the 
total  velocity  gradient 


(5.1) 


(5.2) 
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Superposed  dots  denote  material  time  derivatives,  d/dt.  Both  l}  and  LF  have  components  referred 
to  the  same  intermediate  conhguration  B  introduced  in  Section  2. 

Consider  a  single  crystal  with  one  or  more  glide  systems  labelled  by  index  k.  Let  ^0  denote 
the  line  length  per  unit  reference  volume  of  dislocations  associated  with  system  k,  and  let  and 

denote  the  associated  Burgers  vector  and  unit  normal  to  the  slip  plane  in  the  undeformed  perfect 
crystal.  Restricting  attention  to  slip  systems  with  the  same  Burgers  magnitude,  let  —  bs^^\ 
where  is  a  unit  vector  in  the  direction  of  slip  in  the  perfect  reference  lattice.  Associate  a  uniform 
velocity  with  each  dislocation  population;  immobile  dislocations  have  =  0.  Edge,  screw 
and  mixed  dislocations  and  loops  can  contribute  to  the  same  or  different  densities  depending 
on  their  velocities  and  Burgers  vectors.  Generalisation  of  (3.13)-(3.15)  and  (4.6)-(4.8)  to  multiple 
slip  systems  leads  to  the  usual  rate  kinematics  of  continuum  crystal  plasticity  theory: 

lP  =  ^  0  ^  0  (5.3) 

k  k 


Now  consider  generalised  to  a  volume  element  containing  multiple  dislocations  of  various 
orientations  and  various  edge  and  screw  components.  Motivated  by  developments  in  Sections  3 
and  4,  the  linear  contribution  to  depends  on  the  path  each  dislocation  has  traveled  through  the 
volume  element.  The  magnitude  of  each  dislocation’s  contribution  to  is  proportional  to  its  slipped 

distance  through  the  volume  element  under  consideration.  For  example,  an  edge  dislocation  that  has 
traveled  almost  completely  across  a  cylinder  as  in  Section  3  would  contribute  ^  '^he 

slipped  distance  ^  2R.  It  is  assumed  here  that  the  contributions  of  all  dislocations  comprising  a 
given  family  k  are  equivalent  to  those  of  a  density  located  at  the  centre  of  the  volume  element, 
such  that,  from  generalisation  of  (3.18)  and  (4.9), 

4n  =  E  ®  ^  ® (5.4) 

ut  LyJ  It 


The  factor  is  replaced  by  if  a  hexahedral  volume  element  with  equal  edge  lengths  is 
considered  instead  of  a  cylinder.  Dislocations  with  parallel  Burgers  vectors  of  equal  length  but 
opposite  sign  must  be  assigned  to  distinct  populations  with  parallel  but  oppositely  signed 
directions  (for  example,  24  populations  for  conventional  {111}(110)  slip  in  FCC  crystals).  With 
this  convention,  ^  0  and  ^  0  can  be  enforced.  A  slightly  different  form  for  the  strain 
rate  contribution  from  dislocation  generation  (dislocation  density  rate  distinct  from  slip)  has  been 
proposed  elsewhere  for  shock  compressed  polycrystalline  metals  (29). 

The  non-linear  contribution  to  F^  is  independent  of  the  path  of  dislocation  motion  (the  orientation 
of  the  glide  plane  and  the  slipped  distance),  and  depends  only  on  the  number,  orientation  and  location 
of  dislocations  within  the  element.  For  single  slip  (k  =  1)  of  edge  dislocations,  straightforward 
generalisation  of  (3.29)  leads  to 


-Ulin  ^nonlin  =  " [ln(7r C )  +  l]^(s  ®  S  +  /M  0  /m). 


Ab^ 


(5.5) 


where  A  is  defined  in  (3.22)  and  depends  on  linear  and  non-linear  elastic  properties.  A  similar 
approach  with  (4.14)  and  (4.15)  could  be  used  for  a  single  system  of  screw  dislocations.  An 
analytical  non-linear  solution  to  (3.21)  is  not  available  for  dislocation  networks  often  formed  under 
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activity  of  multiple  interacting  slip  systems;  non-linear  solutions  for  even  isolated  dislocations  are 
scarce  and  approximate  as  mentioned  already.  Motivated  by  forms  of  for  isolated  edge 

and  screw  dislocations  in  Sections  3  and  4,  as  a  further  approximation  it  could  be  assumed  that 
each  dislocation  line  contributes  a  radial  deformation  (a  radial  expansion  or  dilatation  for  material 
properties  considered  herein).  For  many  dislocation  lines  of  various  orientations,  the  net  contribution 
would  then  be  approximated  as  a  pure  volume  change: 

^Llin  ®  (5.6) 

where  ^  is  a  dimensionless  scalar  function  of  dislocation  densities.  Letting  ^ 

<p|f=0  =  0,  0.001 -0.1,  (5.7) 

where  the  estimated  magnitude  at  maximum  dislocation  density  is  motivated  by  Tables  1  and  2  and 
therefore  applies  for  a  representative  crystalline  material  with  properties  in  Table  A.  1 .  Then 

4onlin  ^  /^Lnlin  =  =  E  ^^.8) 

A  dilatational  contribution  to  total  deformation  has  been  confirmed  experimentally  for  several  metals 
(26).  Highly  structured  anisotropic  dislocation  arrangements,  such  as  those  occurring  in  dislocation 
cell  walls,  would  presumably  also  give  a  deviatoric  contribution  to  but  an  analytical  form  for 

such  contrihutions  remains  undetermined.  Combining  (5.4)  and  (5.8)  and  adding  a  core  contribution 
as  in  (3.17), 

=  -^lin  +  ^nonlin  +  ^core-  (5-9) 

where  is  deviatoric  (traceless)  from  (5.4),  is  estimated  as  spherical  (diagonal)  in  (5.8) 

(30),  and  ^  must  he  determined  from  atomic  theory.  The  particular  form  of  the  latter 
is  unknown  from  the  present  non-linear  elastic  analysis,  but  iSj-gj-g  must  vanish  when  ^  =  0  and 
presumably  depends  on  For  an  isotropic  uniform  distribution  of  screw  dislocations, 
could  be  estimated  by  time  differentiation  of  (4.19)  and  (4.20),  but  such  an  approach  would  not  be 
valid  for  edge  dislocations  or  anisotropic  solids,  since  atomic  calculations  (22)  show  that  the  core 
contribution  is  not  spherical. 

The  total  velocity  gradient  of  (5.1)  can  now  be  written 

L  =  fV-‘  [-«■>  +  s<»  ®  *(«  +  +  Z.L, 

(5.10) 

where  quantities  pushed  forward  by  the  total  (recoverable  +  residual)  lattice  deformation  are 
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Elastic  part  ^  is  the  remainder  if  L  is  imposed,  and  is  associated  with  rigid  spin  and  the  rate 

of  stress  applied  to  the  element’s  boundaries. 


5.2  Energy  in  constitutive  models  of  crystals  with  defects 

Let  €  be  a  frame  indifferent  measure  of  finite  elastic  strain,  for  example, 

e  =  1(FE  VE  _  1)^  ^  1  (^E^E  _  (5  14) 


Let  f  be  a  dimensionless  internal  state  variable  accounting  for  effects  of  dislocations  on  energy 
density,  for  example. 


I 


(5.15) 


Total  strain  energy  (free  energy  omitting  temperature  dependence)  per  unit  volume  in  configuration 
Bq  for  an  element  of  crystalline  material  with  dislocations  is  posited  as 


'T  =  'l'(e,  t)  =  W(e)  +  ((I)  +  ^). 


(5.16) 


Letting  ()|o  denote  a  quantity  evaluated  at  the  reference  state  wherein  e  —  0  and  ^  =  0,  a  series 
expansion  is 
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(5.17) 

(5.18) 

(5.19) 


where  constant  and  linear  terms  in  strain  or  defect  density  are  eliminated  so  that  conjugate 
thermodynamic  forces  vanish  in  the  reference  state.  Coefficients  of  strain  in  elastic  energy  W  are 
second-,  third-,  . .  .order  elastic  constants  for  a  perfect  crystal  (see  Appendix  A).  Lunction  i  accounts 
for  self-energy  of  dislocations  within  the  externally  unloaded  (self-equilibrated)  volume  element, 
including  local  elastic  self-energy,  core  energy  and  interaction  energies  (31)  among  dislocations 
within  the  element  (longer  range  or  inter-element  dislocation  interactions  are  measured  by  e). 
Lunction  f  accounts  for  coupling  between  applied  loading  and  defects,  and  while  usually  omitted  in 
constitutive  models,  should  not  be  arbitrarily  dismissed  when  defect  density  approaches  its  theoretical 
maximum.  Lor  example,  the  third  term  on  the  right  in  (5.19)  accounts  for  influence  of  defects  on 
tangent  elastic  coefficients  (32). 

Term  F^  is  not  included  explicitly  in  (5.16)  since  it  is  assumed  that  contributions  from  defects 
are  fully  represented  by  f .  However,  §  =  and  F^  or  its  spatial  gradient  are  not  necessarily 
independent.  Let  C,  —  where  f  and  f  ^  are  the  tensor  density  of  geometrically  necessary 

dislocations  and  scalar  density  of  statistically  stored  dislocations,  the  latter  having  no  net  Burgers 
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vector.  In  the  theory  of  continuously  distributed  dislocations,  contained  in  a  local  volume  element 
is  related  to  that  element’s  plastic  curvature  (15,  25,  33): 


def  ^  ^p^^  (5.20) 

where  —  detF^  =  1  and  V  is  the  material  gradient  in  the  present  context,  ejjK  are  permutation 
symbols,  and  the  algebraic  sign  convention  and  order  of  indices  in  f  varies  in  the  literature  (1,  2, 
25,  33).  It  can  be  shown  that  (34) 


rP —  1  77'P  ^  O  I7P  rE  tI  —  1  I7E —  1  o  /  TT'I —  1  77'E —  1  \ 

J  Fal^UKdjFfiK  =  J  J  F^g  Fg.  eijkdjiFp^  F^^  ). 


Though  not  necessary,  if  ^  obeys  integrability  conditions  d^F^-  ^  ^  ’  '^hen 


>G  tE  tI 1  pE— 1  p-E— 1  0.17I— 1 


which  implies  f  =  ^((^^,  F^,  F^,  VF^)  with  V  here  the  spatial  gradient,  and  thus 

'h  =  4/(6, 1)^4'  =  'h(F^,  F\  VF\  C^).  (5.23) 

Note  that  dkF^^  =  ~Fly} d^F^^.  Global  integrability  of  F^“^  suggests  that  stress  acting  on 

the  boundary  of  each  volume  element  can  be  relaxed  (in  conjuction  with  local  rotation  field  R^~^  ) 
so  that  B  is  a  globally  compatible  conhguration,  presuming  that  the  crystal  is  simply  connected 
in  configuration  B.  Recall  that  in  the  theory  of  continuously  distributed  dislocations,  dislocation 
core  volumes  and  displacement  jumps  across  internal  slipped  surfaces  are  not  resolved  explicitly, 
so  simple  connectivity  of  a  deformed  body  (B)  containing  dislocation  densities  becomes  possible 
(25).  Statistically  stored  dislocations  exert  no  long-range  residual  elastic  helds  associated  with 
incompatibility,  but  still  contribute  to  (5.23)  through  short-range  and  core  energies. 

Typical  gradient-type  crystal  plasticity  models  (2, 15, 35)  invoke  the  usual  two-term  decomposition 
F  =  F^F^  of  (1.1),  and  the  presence  of  dislocations  within  a  local  volume  element  is  indicated  by 
V  X  F^  ^  0  as  in  (5.20).  However,  this  description  alone  does  not  account  for  statistically  stored 
dislocations,  which  may  exist  regardless  of  whether  or  not  the  curl  of  the  plastic  deformation  gradient 
vanishes;  a  non-vanishing  curl  only  indicates  the  presence  of  geometrically  necessary  dislocations. 
In  pure  ductile  metals  under  homogeneous  deformation,  statistically  stored  dislocations  would  be 
expected  to  contribute  most  strongly  to  strain  hardening  and  stored  energy.  Furthermore,  since 
detF^  =  1  in  slip-based  crystal  plasticity  theory,  local  residual  dilatation  cannot  be  captured  by 
such  models  invoking  (1.1).  In  contrast,  and  though  perhaps  more  complicated,  the  present  theory 
invoking  (1.2)  can  account  for  residual  dilatation  (as  well  as  residual  lattice  shape  change)  via  the  F^ 
term,  which  in  turn  is  affected  by  the  total  dislocation  density,  including  both  geometrically  necessary 
and  statistically  stored  dislocations. 


5.3  Constitutive  relations  and  dissipation  inequality 

Under  the  present  idealisation  of  isothermal  conditions,  the  local  balance  of  energy  and  entropy 
inequality  are  (35) 

U  =  t:L,  (5.24) 
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where  U  is  internal  energy  per  unit  reference  volume,  x  —  Ja  is,  macroscopic  symmetric  Kirchhoff 
stress  with  a  Cauchy  stress  and  J  —  detF  =  —  J^J^,  and  'I'  represents  free  energy  density. 

Using  (5.10)  and  (5.16), 

[/  =  T  :  [L^  +  F^F^{LF  + 


=  (F- 


e-ItFE-T)  -.k  +  bY  +-^y 

^  2v7r  ^ 


(5.25) 


is  total  external  stress  power,  and 


a'h  O'h  . 

'h  =  -  :  €  H - 1, 

d€  at  ^ 


(5.26) 


is  the  free  energy  rate,  where  resolved  shear  stress  and  reference  pressure  (energy  per  unit  reference 
volume)  are 

^(k)  ^  .  ^^ik)  ^  \  F  —-Xkk  =  Jp- 


Substituting  into  the  second  of  (5.24)  and  considering  admissible  processes, 

p  pa'h  pT  def  34' 

t(FE,  t)  =  fE— fET,  h(€,  t)  =  — ; 

d€  at 


b  Y  ^  -3fY  +  r  :  L, 

k  k  k  ^ 


cl 

core 


(5.27) 

(5.28) 

/it  ^  0. 

(5.29) 


Equations  (5.28)  are  the  hyperelastic  stress-strain  law  and  definition  of  the  thermodynamic  conjugate 
force  to  internal  state  variable  t.  If  (5.15)  applies,  then  dissipation  inequality  (5.29)  becomes 

k  k 

(5.30) 


\2Jn 


{k)^^{k)y\l2  , 


3F 


dcp 

Yk) 


-hY'^ 

2 


^0. 


Consider  each  term  in  (5.30).  The  first  term  is  dissipation  from  slip  of  mobile  dislocations.  The  second 
is  stress  power  associated  with  deformation  of  the  volume  element  from  core  effects,  noting  that  only 
the  symmetric  part  of  contributes.  Remaining  terms  in  square  brackets  are  associated  directly 
with  dislocation  density  rates  The  first  term  is  due  to  shearing  associated  with  the  rate  of  F^,  in 
the  simple  case  of  an  isolated  dislocation  resulting  from  the  linear  elastic  contribution  to  -  The 
second  is  due  to  dilatation  associated  with  the  rate  of  F^  for  an  isolated  dislocation  resulting  from  the 
non-linear  elastic  contribution  to  The  b^h  term  represents  free  energy  storage  associated  with 

microscopic  strain  energy,  core  energy  and  interaction  energies  of  dislocations  within  the  volume 
element.  Dissipation  from  slip  is  positive  when  resolved  shear  stress  and  dislocation  velocity 
share  the  same  sign.  Remaining  terms  can  contribute  positively  or  negatively  to  dissipation/storage 
of  energy  because  algebraic  signs  of  and  its  multipliers  are  unrestricted;  dislocation  rates  may 
be  positive  or  negative  as  defects  enter  or  exit  the  volume,  or  are  generated  or  annihilated. 
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6.  Examples 

6. 1  Simple  shear 

Consider  a  single  crystalline  cylindrical  volume  element  of  initial  volume  V  =  itR^L  with  L  parallel 
to  unit  vector  e^,.  Let  this  element  be  subject  to  homogeneous  simple  shear  deformation  of  the  form 

F  —  =  1  +  Hci  ®  62  =  1  +  Ss  ®  (6-1) 

Applied  shear  deformation  is  S,  and  the  crystal  is  assumed  perfectly  oriented  for  single  slip,  with 
s  =  6i  and  m  —  €2-  For  this  2D  problem,  dislocation  density  C,  consists  of  a  single  family  of  edge 
dislocations  with  b  —  bs.  Since  s  ■  m  =  0,  cumulative  plastic  deformation  is  of  the  form  in  (3.13): 

F^  =  exp(ys  ®  m)  =  1  +  ys  ®  wi,  (6.2) 

with  cumulative  slip  y  to  be  determined  in  the  analysis  as  a  function  of  S.  A  rate  independent  yield 
criterion  is  used  here  since  S  is  increased  monotonically  at  a  fixed  rate: 

r  T  :  (s  ®  m)  <  To  +  aGh^,  (6.3) 

with  r  the  resolved  Kirchhoff  stress,  tq  the  initial  yield  strength,  and  a  a  constant  in  the  usual  Taylor 
hardening  model  (36).  Linear  hardening  with  modulus  H  is  assumed  such  that  Hy  —  ah^,  leading 
to  the  following  dimensionless  yield  criterion  and  evolution  equation  for  dislocation  density: 

r/G  <  ro/G  +  Hy-  C  =  [Hy/{ab)f.  (6.4) 


Residual  lattice  deformation  is,  from  (3.30),  the  second  of  (6.4),  and  assuming  null  core  traction, 

r2„2 


T  Hy 

F\y)  =  1  +  ( 

yjna 


AH'^y'^  jxH^y^ 

)m - In - r — (s  ®  s  +  m  ®  m). 

4 


(6.5) 


Under  isothermal  conditions,  and  omitting  coupling  between  internal  state  variable  |  =  b^  and 
elastic  Green  strain  €  —  —  1),  (5.16)  and  (5.17)  can  be  combined  to  give  free  energy 

density  'I'  of  the  form  below,  to  order  three  in  elastic  strain: 

'F  =  ~^oifiyS€afi^yS  4“  -^^afiy&€(t>^afi^y6^i<p  4“  t(^)>  (6.6) 

where  second-  and  third-order  elastic  constants  CapyS  and  CapyS(;<p  are  described  in  Appendix  A, 
and  stored  energy  of  cold  work  l  is  inconsequential  in  the  stress-strain  response  in  this  example. 
Kirchhoff  stress  is,  from  (5.28), 


r(fE)  =  E^ - =  (6.7) 

9€ 

with  a  the  Cauchy  stress,  noting  J  —  detF  =  —  1  from  (6.1)  and  (6.2). 

Solution  of  the  problem  proceeds  as  follows.  Deformation  3  is  updated  incrementally.  For  each 
increment,  use  of  (6.1),  (6.2),  (6.5),  and  (6.7)  enables  iterative  solution  of  yield  condition  (6.4) 
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Fig.  5  Shear  stress  (a)  and  Cauchy  pressure  (b)  predicted  using  three-term  (F  =  F^F^F^)  and  two-term 
(F  =  F'^F^)  models  for  single  crystal  in  simple  shear  and  oriented  for  single  slip,  with  hardening  modulus  H 


for  slip  y,  with  inequality  holding  during  elastic  loading  (y  =  0)  and  equality  holding  upon  yield 
(y  >  0).  With  F(3)  and  y  known,  all  other  quantities  of  interest  (cr,  etc.)  can  be  computed. 
Representative  material  properties  are  taken  from  Table  A.  1  of  Appendix  A,  leading  to  A  =  as  in 
Section  3.  Also,  tq  =  otG,  a  =  and  two  values  of  H  and  1)  are  considered. 

Shear  stress  cti2  and  pressure  p  —  —  jtrcr  are  shown  in  Figs  5(a)  and  5(b).  Results  labelled 
‘F  =  F^F^’  are  obtained  by  setting  F^  =  1  instead  of  (6.5),  with  all  other  model  features  the 
same  as  for  the  three-term  theory  with  F  =  F^F^F^  described  above.  Shear  stress  increases  in  the 
plastic  regime  with  increasing  linear  hardening  modulus  //,  as  expected.  Shear  stresses  are  lower 
when  the  three-term  model  is  used  because  F^  accommodates  some  of  the  applied  shear  strain  that 
would  otherwise  be  accommodated  by  F^.  Pressures  are  higher  when  the  three-term  model  is  used 
because  F^  includes  dilatation  7^  =  1  >  1 .  Although  nearly  the  same  shear  stress  response  could 

be  obtained  by  selecting  different  values  of  H  for  two-  and  three-term  models,  only  the  three-term 
theory  is  able  to  properly  account  for  pressure  rise  resulting  from  dilatation  due  to  dislocations;  such 
dilatation  has  been  observed  in  experiments  (26)  and  atomic  simulations  (21,  28). 

6.2  Uniaxial  compression 

Now  consider  a  single  crystalline  cylindrical  volume  element  of  initial  volume  V  =  rcR^L  with  L 
parallel  to  unit  vector  ei .  Let  this  element  be  subject  to  homogeneous  uniaxial  strain  of  the  form 

F  =  F^F^F^  =  (1  —  3)ei  ®  ei  +  e2®  e2  +  ej,.  (6.8) 

Applied  compression  is  S  ^  0,  and  the  crystal  structure  is  taken  as  FCC,  with  cube  axes  [100], 
[010]  and  [001]  parallel  to  e\,  €2  and  ^3,  so  compression  occurs  along  [100]  in  this  problem. 
This  example  is  representative  of  planar  shock  compression  of  a  cubic  single  crystal  (37),  though 
here  temperature/entropy  effects  are  omitted  for  simplicity.  From  symmetry  (37),  8  of  12  possible 
[111}{110)  slip  systems  are  active  upon  yielding  (k  =  1,2,...  8),  each  with  the  same  cumulative 
slip  y  =  y^*^  =  f  y^^Mt.  For  monotonic  loading,  integration  of  (5.3)  results  in  plastic  deformation 
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of  the  following  form,  to  order  three  in  y  ■ 


(6.9) 


with  cumulative  slip  y  to  be  determined  in  the  analysis  as  a  function  of  3.  Equations  similar  to  (6.3) 
and  (6.4)  apply: 

=  r  :  (s^^^  0  /w^^^)  <  tq  +  aGb^  =  tq  +  GHy,  ^  =  [Hy/(ab)]'^.  (6.10) 

Here  ^  again  is  assumed  to  consist  entirely  of  straight  edge  dislocations  with 

=  bs^^\  The  linear  elastic  contribution  to  is  given  by  analogue  of  (3.18)  and  (5.4);  the 
non-linear  elastic  contribution  a  volume  change  of  the  form  in  (3.22),  (3.23),  (5.6),  and  (5.8);  core 
traction  is  again  omitted  for  simplicity.  These  assertions,  with  the  second  of  (6.10),  result  in 


V  OJL  LI  j 

k 


(6.11) 


The  non-linear  hyperelastic  response  is  again  dictated  by  (6.6)  and  (6.7).  The  same  numerical  solution 
procedure  and  material  properties  used  in  Section  6.1  apply  here. 

Axial  true  stress  — crn  and  dislocation  density  ^  are  shown  in  Figs  6(a)  and  6(b),  the  latter 
normalised  by  ^max  =  Results  labelled  "F  —  F^F^'  are  obtained  by  setting  F^  —  \  instead  of 
(6.11),  with  all  other  model  features  the  same  as  for  the  three-term  theory  with  F  —  F^F^F^ .  Axial 
stresses  [Fig.  6(a)]  are  identical  for  all  cases  until  yield  is  attained  at  3  =  1  —  7  Rs  0.017,  which 
would  correspond  to  the  Hugoniot  elastic  limit  in  a  planar  shock  experiment.  Stresses  increase  in 
the  plastic  regime  with  increasing  linear  hardening  modulus  H,  as  expected.  Differences  between 
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Fig.  6  Axial  stress  (a)  and  dislocation  density  (b)  predicted  using  three-term  (F  =  F^F^F^)  and  two-term 
(F  =  F^F^)  models  for  FCC  single  crystal  compressed  along  [100],  with  hardening  modulus  H 
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two-  and  three-term  models  are  negligible  for  the  lower  hardening  modulus  //  =  5,  but  not  for 
higher  hardening  with  H  —  Dislocation  densities  [Fig.  6(b)]  predicted  by  the  two-term  theory 
increasingly  exceed  those  predicted  by  the  three-term  theory  as  y  or  //  is  increased.  At  S  ^  0.2, 
dislocation  densities  predicted  using  //  =  j  are  of  comparable  magnitude  to  maximum  densities 
observed  during  shock  compression  experiments  on  Cu  and  Cu-Al  alloys  (38). 

6.3  Discussion 

The  presentation  in  Sections  6.1  and  6.2  demonstrates  how  a  constitutive  model  incorporating  the 
three-term  decomposition  (1.2)  can  be  implemented,  in  a  straightforward  way,  for  crystal  volumes 
undergoing  homogeneous  deformation.  Results  have  indicated,  for  representative  material  proper¬ 
ties,  the  applied  strain  required  for  the  term  to  significantly  affect  stresses.  Differences  between  the 
present  theory  and  the  usual  two-term  description  become  most  apparent  at  large  strain  in  crystals  with 
strong  hardening  associated  with  defect  accumulation.  Validity  of  the  theory  for  addressing  residual 
lattice  deformation  has  been  established,  within  an  order  of  magnitude,  by  comparison  with  results 
of  atomic  simulation  in  Table  1 .  Comparison  of  theoretically  predicted  dilatation  with  experiments 
reporting  energy  of  cold  work  and  residual  volume  change  has  been  presented  elsewhere  (17,  26), 
wherein  again  order-of-magnitude  accuracy  has  been  confirmed.  As  suggested  by  derivations  and 
calcuations  in  Sections  3  and  4  for  edge  and  screw  dislocations,  the  magnitude  of  residual  lattice 
contribution  from  defect(s)  within  a  local  volume,  F^,  can  exceed  several  percent  as  the  theoretical 
maximum  dislocation  density  is  approached,  but  is  small  for  dilute  dislocation  concentrations.  Here 
‘magnitude  of’  refers  to  the  largest  absolute  value  of  any  stretch  component.  Bulk  dislocation 
densities  of  up  to  2-10  per  cent  of  the  theoretical  maximum  have  been  observed  in  atomic  and 
discrete  dynamic  simulations  of  shock  loading  (39),  which  would  impart  substantial  contributions 
to  F^  with  a  corresponding  effect  on  observed  pressure-density  behaviour  in  shock  compression 
experiments  (30).  Even  if  the  bulk  dislocation  density  is  small  when  averaged  over  a  large  volume  of 
crystal,  the  local  dislocation  density  at  volume  elements  comprising  dislocation  cells  and  cell  walls 
during  latter  stages  of  work  hardening  may  be  quite  large  (for  example,  approaching  the  theoretical 
maximum),  so  the  magnitude  of  F^  would  be  substantial  for  those  elements. 

More  extensive  algorithms  and  numerical  simulations  are  needed  for  applying  the  theory  to 
problems  involving  heterogeneous  deformation;  such  developments  are  reserved  for  future  work. 
For  example,  local  regions  with  dislocation  densities  on  the  order  of  10^^  m~^  have  been  reported  for 
a  nickel  crystal  subjected  to  wedge  indentation  (40),  corresponding  to  ^ /(i’max  ^  2  x  10“^.  Referring 
to  results  in  Table  1,  F^  would  be  expected  to  be  non-trivial  in  such  regions. 

By  definition,  ^  1  as  the  dislocation  content  in  a  volume  element  goes  to  zero  0).  If,  as 
in  the  model  of  Sections  6. 1  and  6.2,  plastic  strain  hardening  results  from  accumulated  dislocations, 
then  F^  ^  1  for  an  elastic-perfectly  plastic  crystal  through  which  dislocations  pass  but  do  not 
accumulate. 

Several  (co)authored  previous  works  (7,  8,  9,  17)  proposed  three-term  decompositions  but  did 
not  derive  their  origins  in  the  context  of  elastic  solutions  for  single  dislocations  and  did  not 
develop  constitutive  frameworks  identical  to  that  in  Section  5.  Content  in  the  monograph  (12) 
mostly  summarises  content  of  such  prior  work,  in  Section  9.1.3  leaving  the  physical  definition  of 
the  third  term  ambiguous  in  general  thermomechanical  analysis,  similar  to  (5).  A  different  three- 
term  decomposition  of  the  deformation  gradient,  written  F^  —  F^^F^^F^,  has  also  recently  been 
developed  (14).  In  this  different  approach,  following  early  concepts  of  dislocation  field  theory 
(1,  4),  fS  =  F  is  the  total  shape  deformation  of  a  material  element,  F^  represents  rigid  slip  due 
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to  dislocations  currently  within  the  element  and  those  that  have  passed  completely  through,  is 
lattice  deformation  that  returns  the  element  to  a  compact  state  (residual  elastic  and  core  effects),  and 
is  recoverable  lattice  deformation  associated  with  loading  of  the  external  boundary  of  the  entire 
body.  Generally,  F^  and  F^^  are  incompatible  fields  (they  do  not  obey  integrability  conditions), 
but  is  compatible  and  is  the  spatial  gradient  of  a  vector  field.  Both  the  present  approach 

and  that  in  (14)  derive  a  multiplicative  decomposition  for  single  edge  and  screw  dislocations,  and 
then  generalise  to  multiple  defects.  However,  steps  followed  in  each  approach  differ,  leading  to 
different  interpretations  of  the  third  (middle)  term  in  the  decomposition  associated  with  residual 
lattice  distortion.  This  leads  to  different  numerical  results  in  the  present  work  and  in  (14).  Example 
calculations  for  simple  shear  and  uniaxial  compression  in  Sections  6.1  and  6.2  are  the  first  such 
reported  using  the  current  three-term  description  and  have  not  been  presented  elsewhere.  Appendices 
herein  contain  non-linear  elastic  derivations  and  treatment  of  anisotropy  not  given  in  (14). 

7.  Conclusions 

A  three-term  multiplicative  decomposition  for  an  element  of  single  crystal  containing  dislocations 
has  been  derived.  In  addition  to  usual  terms  associated  with  elastic  deformation  and  cumulative 
dislocation  slip  through  the  element,  a  third  term  has  been  incorporated  that  accounts  for  defects 
within.  This  third  term  is  further  divided  into  contributions  from  the  incompletely  slipped  surface 
(a  ‘linear’  term),  elasticity,  and  core  effects.  Expressions  have  been  derived,  and  magnitudes 
estimated,  for  contributions  to  this  third  term  for  cylindrical  elements  containing  a  single  edge  or 
screw  dislocation  line.  For  dislocation  densities  approaching  the  theoretical  limit,  such  magnitudes 
may  be  comparable  to  applied  elastic  strains  needed  to  initiate  slip.  Upon  generalisation  of  the 
theory  to  crystals  with  multiple  defects,  appropriate  general  forms  of  free  energy  functions  have 
been  suggested  and  analysed  in  the  context  of  the  three-term  decomposition.  The  theory  has  been 
implemented  in  constitutive  model  calculations  for  simple  shear  and  uniaxial  strain  compression, 
with  differences  in  results  from  those  of  usual  two-term  crystal  plasticity  apparent  for  large  strains 
and  strong  plastic  hardening.  Primary  advantages  of  the  proposed  three-term  description  are  (i)  its 
precise  delineation  of  contributions  from  dislocation  motion  and  dislocation  generation  to  irreversible 
deformation,  and  (ii)  its  ability  to  quantify  local  residual  elastic  strains,  including  volume  changes. 
Neither  phenomenon  can  be  captured  by  two-term  crystal  plasticity  wherein  the  rate  of  plastic 
deformation  results  only  from  slip  rates  described  by  the  usual  Orowan  equation. 
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APPENDIX  A 

Non-linear  elastic  coefficients 

Let  Ua(Xp)  denote  elastic  displacement  and  Xfj  Lagrangian  coordinates;  Green  elastic  strain  is 


/  dug 

[dXp 


dujj 


dliy  Blly 


1 

Egp  -  - 


(A.1) 
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Strain  energy  density  W  written  as  a  Taylor  polynomial  is 

^  T  '^^ajiiyS€(p^alii^y8^e<j>  +  *  *  ‘  •  (A. 2) 

For  homogeneous  isotropic  materials,  elastic  constant  tensors  are  of  the  form 

^ctfiyS  —  K^apSys)  +  G(SaySps  +  SasSpy),  (A. 3) 

^afiyS€<p  — 

+  V2[ScipiSy(;Ss(l)  +  Syt/iSsc)  +  ^y8i&ae&P4,  +  ^acji&pa)  +  S^^(SaySps  +  ^asSpy)] 

+  V^l^ayiSpeSs^  +  Sp^S^f:)  +  SpsiSafSy^  +  &a<l>^ye) 

+  ^aSi^PeSy^  +  Sp^Sya)  +  SpyiSafSs^  +  Sa^Ss^)].  (A.4) 

Let  C/7  and  C[jk  denote  second-  and  third-order  constants  in  Voigt  notation,  where  indices  1,  2,  ...  6.  Second- 
order  elastic  constants  obey  the  familiar  relations 

L  =  Ci2,  G  =  C44  = -(Cii  -  C12),  K  =  X+-G,  V  =  2(k  +  G)’ 

where  G,  K  and  v  are  shear  modulus,  bulk  modulus  and  Poisson’s  ratio  in  the  reference  state.  Third-order  elastic 
constants  obey 

Vl=Ci23,  V2  =  Ci44  = -(C112  -  C123),  V3  =  C456  =  g(Clll  -  3Cii2 -F  2C123).  (A.6) 

Let  ^(p)  and  G(p)  denote  incremental  bulk  and  shear  moduli  of  an  isotropic  elastic  solid  subjected  to  dilatation 
induced  by  Cauchy  pressure  p.  Then  (15,  41) 

K'  "  (d^/dp)|p=o  =  -[vi  +  2v2  -F  (8/9)y3]/^:,  (A.7) 

G'  =  (dG/dp)|p=o  =  -[L  -F  G  -F  V2  -F  (4/3)y3]/^:.  (A.8) 

Ultrasonic  measurements  and  various  high  pressure  experiments  suggest  that  the  pressure  derivative  of  the  bulk 
modulus  in  the  reference  state  obeys  2  <  K'  <1  for  most  pure  elemental  crystalline  solids,  with  3  <  K'  <  6 
typical  (42).  Experimental  data  also  suggest  that  the  pressure  derivative  of  the  thermodynamic  shear  modulus 
in  the  reference  state  obeys  0  <  G'  <  3  for  most  pure  elements  (41).  Representative  values  used  in  this  article 

are  tabulated  in  Table  A.l.  Note  X  —  G  O  v  —  whereas  for  crystalline  metals  and  minerals  typically 

0.05  <  y  <  0.45. 

Table  A.l  Representative  elastic  constants  and  Burgers  magnitude  b. 


XjG 

K' 

G' 

V2/G 
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11 
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APPENDIX  B 

Radial  displacement  for  screw  dislocation  in  non-linear  isotropic  cylinder 

Expressions  for  cj,  C2  and  C3  of  (4.12)  are  listed  or  newly  derived  here.  Considered  only  is  the  case  when 
pressure  p  vanishes  on  the  dislocation  core  boundary  at  p  =  r.  Constant  C3  is  (15) 


C3  = 


M(G+ V3) 

87r2(X  +  2G)  ’ 


(B.l) 


with  A  the  Lame  modulus  and  V3  a  third-order  elastic  constant.  Constants  ci  and  C2  are  determined  from 
vanishing  traction  conditions  at  p  =  r  and  p  —  R,  requiring  simultaneous  solution  of 


2G  ,  ,  fc2(G+v3)  ^  Inr  ,  ^2^  +  ^2)  _  ^ 

- “I - - ^^3  — o — I - : — —  —  0, 


2G  ,  ,  b^(G -I-V3)  ^  InR  ,  fc2G  +  v2) 

- T^l  +  G)c2  H - - Gc3  — 3 — - 3 — T — 

r2  Stt^R^  Sn^R^ 


=  0, 


(B.2) 

(B.3) 


where  V2  is  another  third-order  elastic  constant.  After  some  algebra. 


Cl  = 


C2  = 


167r2 


1-1- 


A  -f  V2  -f  V3  (G  -I-  ^3)  (1  -  r2//?2)in  ^  _  (r2/R2^ln(R/r) 


(A  -f  2G) 


1  -  r^/R^ 


b^ 

r  G(G-fV3)  1 

Inr-  (r2/R2)in(b?/r) 

167r2r2 

.(A -b  G)(A -b  2G). 

1  -  7-2 /r2 

(B.4) 


(B.5) 


While  generic  forms  of  (B.2)  and  (B.3)  are  listed  by  Teodosiu  (15),  particular  solution  (B.4)  and  (B.5)  is  newly 
derived  in  the  present  work.  Choosing  typical  values  from  Table  A.  1, 


Cl 


C2  ' 


807r2 


(r^/R^)ln(R/r)  -  (1  -  r^/R^)lnr 

1  -  r2/f?2 


160jr2r2 


(r^/R^)ln(R/r) 
1  -  r^/R^ 


-5 


C3 


40jr2 


(B.6) 


(B.7) 


Radial  displacement  (4.12)  at  p  =  R  becomes 


Cl  InR 

niR)  =  —  +  C2f^-I-C3  — , 

K  K 


(B.8) 


where  results  in  Fig.  3  and  Table  2  of  Section  4  correspond  to  r  —  b  of  Table  A.l. 


APPENDIX  C 


Elastic  anisotropy 

Closed-form  analytical  elastic  solutions  for  dislocations  in  anisotropic  non-linear  elastic  solids  do  not  appear  to 
exist,  though  numerical  integration  of  lengthy  and  complex  approximations  available  in  the  literature  (19)  may 
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be  possible.  Analytical  solutions  can  be  obtained,  however,  for  volume  change  AIZ/V  resulting  from  (3.21) 
in  crystals  with  cubic  symmetry,  effectively  replacing  the  first  equality  in  (3.22)  or  (4.10).  In  what  follows, 
effects  of  anisotropy  on  dilatation  are  quantified  for  several  kinds  of  straight  dislocations;  presumably,  effects  of 
anisotropy  on  deviatoric  parts  of  are  of  comparable  magnitude,  though  this  cannot  be  proven  at  present. 

For  a  density  f  =  L/F  of  edge,  screw,  or  mixed  dislocations,  the  residual  volume  change  in  a  cubic  crystal 
can  be  written  as  (16, 17) 


AF 


Uk’  -\)e  + 

K 


1 

G 


(C.l) 


where  any  three  independent  second-order  elastic  constants  at  the  reference  state  are  related  by 

i(Cn+2Ci2),  G=^(C„-Ci2),  Z  =  C44-G,  (C.2) 

and  Z'  is  the  pressure  derivative  of  tangent  stiffness  corresponding  to  Z  at  the  reference  state.  Dilatational 
energy  per  unit  dislocation  length  e  and  deviatoric  energies  e'  and  e  are  (16),  (17) 


-e=  —  ( 

2Ljy\dXa 
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e  =  — 

L  JV 
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1  ^  dua  Sup  ^  dua  Sup  ^  1  ^  dua 
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dV, 
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(C.4) 


Total  dislocation  energy  per  unit  length  e  can  be  written  (43) 


,  R 

e  =  e-\-  e  -3-  e  =  k  —  In  — , 
4:rr  r 


(C.5) 


where  k  is  the  energy  factor  that  depends  on  elastic  constants  and  dislocation  geometry.  Define 


/If  i/e,  /IfeVe,  ,X+/+,X=1. 


(C.6) 


Furthermore,  let  /cjso./iso  corresponding  quantities  in  an  isotropic  elastic  solid  containing  a 

dislocation  having  the  same  geometry  and  elastic  properties  (G,  K,  G’ ,  K’),  noting  that  (Z,  Z’ ,  e,f)  all  vanish 
for  the  isotropic  case.  Then  the  ratio  of  residual  elastic  volume  change  in  a  cubic  crystal  to  that  in  such  an 
isotropic  solid  can  be  written 

AF  f{K' -l)/K+f(G' -G/K)/G+f(Z' -Z/K)/Z  K  k 

A,o(.K'  -  1)/K  +ff^^(G'  -  G/K)/G  ^iso  ^iso 

Iff  ^/iso  ^rid/  ^  ^  are  assumed,  then  2  are  reasonable  bounds  on  M  when  considering  the  non-linear 

elastic  properties  reported  for  10  cubic  metals  in  (26). 
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Fig.  C.l  Energy  factor  for  cubic  crystal  normalised  by  energy  factor  in  isotropic  solid 


Considered  now  are  mixed  dislocations  with  tangent  lines  parallel  to  [101].  For  a  ^(llO)]!!!)  dislocation 
(typical  in  FCC  crystals)  and  for  a  ^  (111){110)  dislocation  (typical  in  BCC),  respectively  (43), 


K  =  -{k\  +  2k2  +  K-i)  (FCC), 


/Cl  —  (  C  + 


2v 


K2  =  Kl 


1  -  2v 
2-2u 
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(BCC); 

(C.8) 

G, 

(C.9) 

2v  /  1  \1 

1/2 

^  +z 

(C.IO) 

2v  \  1  —  2v  / 

Ratio  (c/zciso  depends  only  on  dislocation  type,  v,  and  z,  noting  that  z  —  I  for  isotropic  elasticity.  Values  of 
K/Kiso  are  shown  in  Fig.  C.l  over  the  range  0  v  ^  0.49  for  anisotropy  ratios  z  —  ^  and  z  —  4,  properties 
which  hound  many  cubic  metals  (for  example,  z  sa  1.0  for  W  (44),  z  2.4  for  Fe  (28),  and  z  ^  3.2  for  Cu 
(45)),  the  alkali  metals  being  a  notable  exception  (45).  From  Fig.  C.l,  j  <  iz/Kiso  ^  3.  Combining  this  result 
with  previously  stated  hounds  on  M  yields,  from  (C.7),  the  estimate  g  <  AV/AViso  <  6.  Relevance  of  the 
isotropic  approximation  used  in  Sections  3  and  4  to  determine  non-linear  elastic  contributions  to  can  be 
assessed  for  a  particular  cubic  crystal  by  substitution  of  its  elastic  properties  into  the  above  equations,  or  into 
analogous  equations  that  can  be  derived  using  the  same  procedure  with  anisotropic  energy  factors  available  for 
other  edge,  screw,  or  mixed  dislocation  types  (43). 
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